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ABSTRACT 

We derive the Feynman rules for the graviton in the presence of a flat Robertson- Walker 
background and give explicit expressions for the propagator in the physically interesting cases 
of inflation, radiation domination, and matter domination. The aforementioned background is 
generated by a scalar field source which should be taken to be dynamical. As an elementary 
application, we compute the corrections to the Newtonian gravitational force in the present matter 
dominated era and conclude - as expected - that they are negligible except for the largest scales. 



1. Introduction 



Of the non-trivial spacetime backgrounds, those of particular interest in cosmology 
are spatially homogeneous and isotropic. This is a direct consequence of the cosmological 
principle which states that there is no preferred position in the universe. It can be shown 
that the most general background satisfying the above requirements is the Robertson- 
Walker one [1]: 



It is characterized by a function a(t) and a constant k which can be chosen to take on the 
values +1, 0, or —1. The function a(t) - known as the scale factor - is a measure of the 
"radius" of the universe, while the constant k determines the spatial curvature. 

By far the most natural explanation of the observed homogeneity and isotropy in the 
present universe is the assumption of an inflationary era in the evolution of the universe. 
During that era, the physical distance between observers at rest on fixed spatial coordinates 
increases superluminally. Accordingly, the observed universe was once so small that causal 
processes could have established an initial thermal equilibrium across it. This accounts 
for the fact that the cosmic microwave background radiation from different regions of the 
sky is seen to be in thermal equilibrium to within one part in 10 5 [2]. Since any decent 
amount of inflation * redshifts the spatial curvature to insignificance, we can restrict our 
study to the k = geometries: 



Such geometries do not represent solutions of the Einstein equations in vacuum except 
in the case of flat space where the scale factor is equal to one. In all other cases a 
non-trivial stress tensor must be present. Homogeneity and isotropy constrain the stress 
tensor to have only diagonal elements consisting of the density p(t) and pressure p(t): 
By decent we mean at least the 60 e-foldings of inflation needed to explain the observed degree of isotropy. 




(1.1) 



ds 2 = — g?£ 2 + a 2 (t) dx ■ dx . 



(1.2) 
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Tjj,u = [ p(t) p(t) p(t) p(t) ]. This form can arise either from a non-zero temperature matter 
theory or a non-trivial scalar field in zero temperature quantum field theory. To preserve 
the isotropy of space, the matter fields must be scalar since the Lorentz transformation 
properties of all other fields - like fermions or gauge bosons - exhibit preferred directions. 
The determination of the scalar field action which supports an arbitrary spatially flat 
Robertson- Walker background is given in Section 2. 

It is best to set-up the theory on an arbitrary spatial manifold of finite extent - for 
instance, T 3 x 3? - in which case Ax 1 < R. We shall find it very convenient to utilize the 
conformal set of coordinates: 

ds 2 = — dt 2 + a?{t) dx ■ dx = fi 2 (??) (—drf + dx ■ dx^j . (1.3a) 

The relation between the two coordinate systems is given by: 

dt = fl(r]) df] ; a(t) = 0(r/) , (1.36) 

so that a generic power law in co-moving coordinates: 




(1.4a) 



takes the following form in conformal coordinates: 

0(r ? )=m T ^, (1.46) 

with r]{) = t{) (1 — s) _1 . The physically most distinguishable cases consist of the inflating 
universe for s = +oo , the radiation dominated universe for s = \ > the matter dominated 
universe for s = | , and the flat universe for s = 0. We shall also assume throughout that 
Ax R and At] <^ R since then the integral approximation to mode sums which appear 
because we work on T 3 is excellent. * 

In conformal coordinates light reaches the distance Ax = R in time Ar] = R. 
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In order to study gravitational effects not associated with the smallest of scales in the 
class of backgrounds (1.2), we must develop the appropriate perturbative tools - this is 
accomplished in Section 3. As an elementary application of the perturbative results, we 
calculate in Section 4 the response of the gravitational field due to a point source in a matter 
dominated universe. The resulting corrections to the Newtonian long-range gravitational 
force are found to be negligible as expected. Our conclusions comprise Section 5. 

2. Determination of the source background 

The dynamical system under consideration has a background action which consists of 
two parts: 

S = S e + S m . (2.1) 
The gravitational action Sg is obtained from the Einstein Lagrangian: 

^ = Tab R ^ • (2 ' 2o) 

while the matter action <S m is based on a generic scalar field Lagrangian: 

£ m = -\ d V ip g^V^g - v(<p) V=g ■ (2.2b) 

We obtain the following equations of motion for the gravitational field: 

G^ = (8nG)T^ , (2.3a) 

and for the scalar field: 

d»(V=g<r d„ip) = • (2.36) 



The Einstein tensor G^ v is determined by the geometry (1.2): * 

G°° = 3^V ; & = -gV( 2^ + % ) , (2.4) 



a J \ a a? 



Our metric has spacelike signature and R^ vpa = p + ... Greek indices run from to 3, while Latin 
indices run from 1 to 3. Dots as superscripts denote differentiation with respect to the co-moving time t. 
We are also using units where c = % = 1. 



and the covariantly conserved stress tensor T^ y by the matter sector of the theory: 

T ^ = -j= = ( 9 m 9 vP - 9 aP ) d a <p dptp - gT V{fp) , (2.5a) 

T^- v = . (2.56) 
The resulting density p(t) and pressure p(t) take the form: 

T 00 = p=|^ 2 + V($) , (2.6a) 

r3= pg ii = g i3 (\f-v($)) . (2.66) 

In order not to disturb the spatial homogeneity and isotropy of our background geometries, 
we have taken the scalar field to be a function of time only: <p = <j)(t). 

Using (2.4) and (2.6), the gravitational equations (2.3a) become: 

2 1 / a a 2 



By satisfying equations (2.7), we determine the source parameters <p(t) and V(<p) - or, 
equivalently, p(t) and pit) - which support the flat Robertson- Walker background (1.2). 
In conformal coordinates the above equations take the form: * 

2 I/O" 20 /2 



i / n" q' 2 



where the primes over O indicate differentiation with respect to the conformal time rj. It 
remains to be shown that the above choice for the source parameters is consistent with the 
scalar equation of motion (2.3b): 

~ 3a ~ dVU) n 
+ — ^ = . 2.9 

a (30 



x 2 

* Notice that in the case of de Sitter spacetime, s — > +oo, equations (2.7-8) have the proper limit: <f> 
and V{$) -» (8ttG) A" 1 , where A = 3% 2 . 



This is a direct consequence of the conservation equation (2.5b): 

P=--(P + P) ■ (2-10) 
a 

By substituting (2.6) in (2.10) we obtain: 

s(i* a = — (2 ' n) 

which is identical to (2.9). 



3. The Feynman rules 



3.1 The Lagrangian. 



The quantum fields are the graviton h^y(x) and the scalar (f>(x): 



9jJLV — 9[IV + kH^v 



(3.1a) 



+ </> • 



(3.16) 



It turns out that - just like the case of de Sitter spacetime [3] - it is most convenient to 
organize perturbation theory in terms of the "pseudo-graviton" field, ip^y(x), obtained by 
conformally re-scaling the metric: 

g^y = fi 2 g^ v = O 2 (jifju, + Kipni^j ■ (3.2) 

As usual, pseudo-graviton indices are raised and lowered with the Lorentz metric and 
ft 2 = 16ttG. The total Lagrangian is: 

C = -1 R - \ dyip g^V=9 ~ V(ip) ^9 . (3.3) 
By substituting (3.1b) in (3.3) we get: 



C — 7T R 



77=0 



n\ dcp 



-g , (3.4) 



and we can proceed to organize C according to the number of fields present. 



The pure metric part of C: 



£ (0) - ^2 R 



in terms of the re-scaled metric is: 



(3.5) 



'(0) 



K 



n z r + 2 



(-l + g°°) (W' + fi' 2 ) 



-9 > 



(3.6) 



where we used (2.8) and ignored surface terms. After performing many partial integrations, 
(3.6) can be cast in a form identical to that of de Sitter spacetime [3]: 



9/ 



-\^9 pa r V rltpw^m,* ■ (3-7) 



a /o2^ 



All the self-interactions of the graviton in the presence of an arbitrary scale factor Ct(r]) 
can be obtained by expanding expression (3.7). The quadratic part of is: 



.(2) 

'(0) 



ap,p 



+^ ^ n n' , (3.8) 



where ip = ip^. 

The part of C linear in 4> is simpler: 



£(i) = -0 dv<i> 9°"V^g - V^g , 



dV 



K 



< du<P 9 P ° + £' <t> 
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(3.9a) 
(3.96) 



In the last step we used the scalar background equation of motion (2.9) to derive the 
identity: 

<3Ma) 



where we define £(77) as: 



(3.106) 



The following quadratic part emerges: 



= £ ^0 v M ° + \ (£0)' </> 



(3.11) 



We shall also need the part of C which is quadratic in </>: 



d 2 V 



c {2) = -\d^d v <t> g ^^g - ?q^($) ^ V^g , 



lr>6 t-1 / o-4 tf / \/ 



(3.12a) 
(3.126) 



and the associated piece which is quadratic in the quantum fields: 



= _ 1 2 ^ ^ + 1 6 r 1 ( -4 e / y J 



(3.13) 



In order to properly quantize the theory and calculate the various propagators, we must 
fix the gauge. We shall accomplish this by adding the following term to the Lagrangian: 



£ G F = ~Wv F» F v , 



(3.14a) 



where: 



F^ = Q i>^ v - \Sl V> - 20' ^0 + r/ M o £ 



(3.146) 



The resulting gauge-fixed quadratic Lagrangian is: 



(2) (2) J2) M) , 
^GF "^(0) + ^(1) + ^(2) + ^G F 



- &n 



«2 



O + 02 



o ^ + n 



«2 «" 



^ V'mo + V'oo 



O 



a 2 + ^ + 0- 4 £ 2 + o 4 r 1 ( n- V )' 



n 



(3.15) 
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In terms of a kinetic operator we get: 

£ gf = + V" •/> + V-oo n 



a 2 + ^- + o- 4 e 2 + o 4 r 1 (fi- 4 e / )' 



(3.16) 



where we have: 



IT (P T 



2°^ <V °o 



Da 



+ V ^ C } D S + V So" D B ■ 



(3.17) 



Parenthesized indices are symmetrized and a bar above a Lorentz metric or a Kronecker 
delta symbol means that the zero (i.e., rf) component is projected out: 



jr v e- V X Q X V 



The quadratic operators defined in (3.17) are given by: 



D A = ft 



n2 



<9 2 - 



+ 2- 



/2 



(3.18) 



(3.19a) 



(3.196) 



Notice that is the kinetic operator for a massless, minimally coupled scalar in the 
presence of the gravitational background (1.3). 

3.2 The Diagonal Variables. 

The quadratic Lagrangian (3.16) should be brought into diagonal form. There is mixing 
in the first term of (3.16) between ipn and ^oo : 

D^ a i) pa = \ I D A fa} - \{ip u - ifj Q0 ) D A (ipjj - ifj 00 ) - ^00 Da Voo 

- D^ TpOi + ^00 D B V'oo } , (3-20) 
which is removed by transforming to the field variable 

D^ a ^ P a = \{ Idj Da Cij - Idi Da Cjj - Co* D B Co* + Coo D B Coo } 



\ I Cij Da Crs 



\^i{r $s)j ~ <^sj + Coi D B (Or 

+ Coo d b Coo } > 
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-S. 



IT 



(3.21) 



where we defined: 



ipij = dj + Sij Coo 



^oo = Coo 



(3.22) 



There is also mixing between Coo an d <fi- The part of the quadratic Lagrangian involved 



is: 



^mixing = \ COO D B Coo + CoO ft 



n <f> 



+ & + tL. + fl-4 C 2 + Q 4 r 1 ( Q -4 e / y 



(3.23) 



We can always diagonalize this system but, for most scale factors ft (77), the resulting kinetic 
operators will be non-local. In this case the mode functions obey fourth order differential 
equations. However, for a class of scale factors which includes the generic power law (1.4), 
the system can be diagonalized without sacrificing locality or the second order character 
of the mode equations.* In this case, £ m i x i ng takes the form: 



^mixing - I COO ft 



d 2 + 



1 1 



(1- S ) 2 v 2 
q2 , 2s 2 -3s+2 



ft Coo + Coo ft M ^ ft 



1 



(l- s ) 2 ^2 



. 



Let us call the diagonal variables x an d v: 



(3.24) 



Coo = cos 9 x - sin 9 v , 
= sin 6 x + cos 9 v , 



(3.25a) 
(3.256) 



and demand that the part of ^mixing i nv °l vm g X v vanish: 



C xv = 4 X ft { sin(2fl) + cos(26>) } ft u = . 



Of the two available solutions, it is sufficient to consider one of them: 



cos 2 9=j^ s 



sin 2 e = ^- s 



sin 9 cos 9 = — , 



(3.26) 



(3.27) 



* The condition that O(ry) must obey is: — O 2 Q'^j = tan 2 ^^/cic 2 + 0^, where ci,C2 are 

constants and is an angle. 
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since the other merely interchanges the role of x and v. The result is the following diagonal 
form: 

,2 <^ 2 -- 1 



r s - 1 v o 

^mixing ~~ 2 A " 



_|_ 2s 2 — s 



O x + ^ o 



The complete diagonal quadratic Lagrangian is: 
.(2) 



01 + IT^f ^2 



77- 



(3.28) 



^GF ~~ D A (rs \&i(r S s )j ~ &rs + y(oi Co 

+ 1 vD s cV+ 1 xD s aX ? 



where: 



n 



d 2 2 S 2 - S 



d 2 + 
d 2 + 



(1- 


-s) 2 




s 


(1- 




2 


— s 



1 

r? 2 



(3.29) 

(3.30a) 
(3.306) 
(3.30c) 



It becomes transparent that the "^4" graviton modes are six and have purely spatial po- 
larizations, the graviton modes are three and have mixed polarizations, and the "C" 
graviton mode is one with purely temporal polarization. Clearly, the two physical graviton 
polarizations are transverse traceless "^4" modes. 

3.3 The Propagators. 

Manifest spatial translation invariance allows us to write the general linearized solution 



as the following superposition 

d 3 k 



0(77, x) = 



(2^ 



exp[i£ • x] I #(77, fc, A) a^(k, A) + **(t/, fc, A) 4„(fc, A) } , 



d 3 k 



(2tt)< 



exp[ik ■ x] I *(?7, k) a(k) + **(t/, fc) a^(fc) | , 



(3.316) 



where /c = and the index A spans all available polarizations. The free vacuum of the 
theory is defined as the normalizable state that obeys: * 



0) = 



(3.32) 



* Besides the a M „(fc, A), the analogous ghost operators annihilate |0) as well. 
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The mode functions are canonically normalized: 



(77, k, A) **'(t/, fc, A) - **(?7, ife, A) #'(?7, ife, A) = iQ~ 2 VA , 



and are annihilated upon acting the appropriate quadratic operators. In our case: 



(3.33) 



Df *(t/,A;,J) = 



(3.34) 



The pseudo-graviton and scalar propagators are defined by: 



(x;x) = ^o| t|^i/(x) VVxtOe')} 0^ 
iA(x;x ; ) = ^()| r|0(x) |()^ , 



(3.35a) 
(3.356) 



but we must consider the mixed propagator as well: 



(x;x') = (o\t[iP^(x) (j>(x') 



(3.35c) 



Naturally, it is the propagators of the diagonal variables that we shall evaluate at first. 
From (3.29) we conclude that: 



0\T[Cij(x)Crs(x')} |0 

o|T{Coi(x)Cor(^)} |0 

o| r{u(x) |o 

o|t{ x (x) x (^)} |o 



= iA^(x; a/) 2 fy r 6 8 )j - 5ij <5 rs ] , 



= iA 8 B (x;x') [-& 
= iA^(x;x') , 
= iA^(:r;a/) , 



(3.36a) 
(3.366) 
(3.36c) 
(3.36a 1 ) 



where the propagators iA|(x; x') satisfy: 

Df iAj(x;x') =i5^\x-x') 



I = A,B,C 



(3.37) 



We transform back to the original fields with relations (3.22) and (3.25): 



o| r{v>ij(x) iprs(x')} |()^ = zA^(x;x') 2 5 z(r o s)j - S rs +iA s c (x;x') ^ <5 rs 
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+ iA s A (x;x') Sij S rs , 



T^oi(x) ip 0r (x 
T^oo(x) ifj 00 (x 
T^oo(x) i) rs {x 
o| T^ipij(x) 4>{x 
0| T{ifj 00 (x) 0( 
o| t{ 0(x) 0(x 



zA^(x;x ; ) 



ir 



= iA s A (x;x')^- s +iA s c (x;x') T ^- s , 
= iA^(x; x) 5 rs + iA s c (x; x') S rs , 
= - iA\{x- x') ^ Sij + iA s c (x; x') ^ 8^ , 
= -iA A (x;x')^- s +iA s c (x;x')^- s , 
= iA s A (x-x') j^- s +iA s c (x-x') , 



(3.38a) 
(3.386) 
(3.38c) 
(3.38d) 
(3.38e) 
(3.38/) 
(3.38^) 



To assemble expresssions (3.38) in covariant form, we use (3.18) to derive an identity: 



5, {a 8® = sf + Sf 6® + 6° 5° 6 « ^ 



^pv = ip~pn + 2o (ai ° V f )o + V 5 " ^00 , 



0x0 



which allows us to write the pseudo-graviton propagator as: 



(3.39a) 
(3.396) 



nvA s pa 



(x;x) = iA s A (x;x') { ^„T f 
+ iA s c (x;x') 



A 
pa 



+ 



rpC 

pv x pa 



+ iA s B (x; x 1 ) 



pv-Lpa 



t: 



c 



Hv-Lpa 



(3.40) 



where: 



rpA 

pv 1 pa 

rpB 

pv-Lpa 

rpC 

\^v ± pa 



^p(p ^a)y Vpv Vpa 



There is also the mixed propagator which equals: 



pv L 



(x;x')=[-iA s A + iA s c )^- s 



(3.41a) 
(3.416) 
(3.41c) 



(3.42) 
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3.4 The Explicit Form of the Propagators. 



The condition (3.32) on the vacuum and the free field expansions (3.31) imply that 
d s k 



iAj(x; x) 



exp(— ek) exp ik ■ (x — x')j | 9{r] — ?/) ^(77, k, I) ^*(rj' , k, I) 

I = A,B,C . (3.43a) 



(2tt) 3 

+ e( v '- v ) **( v ,k,i)*(rf,k,i) } 



It is elementary to perform the angular integrations: 

1 f +OG r 
iA I (x;x) = — tt— I dk exp(-efc) k sin(fcAx) <^ 9(Arj) ^(77, fc, J) **(?/, fc, I) 
Z7T Aa; Jo ^ 

+ 0(-Atj) **(77, fc, J) fc, /) } , I = A,B,C , (3.436) 

where Ax = ||x — and A77 = 77 — 77'. Notice that we have included the traditional 
ultraviolet convergence factor exp(— ek) which serves as an ultraviolet mode cutoff. 

The actual computation of the propagators is a three-step process. First, we solve 
equation (3.34) to obtain the mode functions ^(77, k, I). Then, we evaluate the propagators 
iAj(x; x') by inserting the respective mode functions in (3.43). In the final step, we simply 
substitute the derived expressions for iAj(x; x') in (3.38g) , (3.40) and (3.42) to arrive at 
the scalar, pseudo-graviton and mixed propagators respectively. 



Step I. Consider, to begin with, the "A" modes equation: 



d*l/(7/,x) = o n 

Now rewrite (3.44) as follows: 

O 2 



-dt + v z + 



1 



2 , 2s 2 -s 



1 



Qf(rj,x) = Q . 



(3.44) 



^ 2 + V 2 -t^-^ 



f(v,x) = 



(3.45) 



and Fourier transform in momentum space: 



d2 rj + rh- d v + k 2 ] /(77,£) = o , 



This can be cast in the form of the Bessel equation: 

/ d 2 Id v 2 \ n 

V dy z y dy y z J 
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(3.46) 



(3.47) 



by substituting: 



y = krj 



g(y) = y w G(y) = f(r I ,k) , 



(3.48) 



so that (3.46) becomes: 



y 



w 



(3.49) 



This is the Bessel equation for w 
functions: 



l-3s 



2(1-*) 

G^ 2 )(y) = 4i ,2) (y) 



The two solutions to (3.49) are Hankel 



V A = 5 



l-3s 



and are complex conjugates of each other: 



A similar analysis is applied to the case of the modes: 



-f) 2 + V 2 -I £ 



In momentum space equation (3.52) becomes: 

f(v,k) = , 

which is again equal - using (3.48) - to a Bessel equation when w = 
r d 2 

111 L«/ 



a 2 + ft^,-ft^ + * 2 



l-3s_ 



+ ( ^ + 2w) ^ + 1+ (" ( ^ + w - 1) -^)^] G(?/) = • (3 - 54) 



<hr ' y ay ' > ir 

The following set of complex conjugate solutions are implied: 



G^\y) = H^\y) 



1+s 



Finally, the "C" modes lead to the following equation: 



B 2 , 2s 1 8 2(1+3) 1 ,2 



/(»7,fe) = • 
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(3.50) 



(3.51) 



n f(n,x) = o . (3.52) 



(3.53) 



(3.55) 



(3.56) 



The resulting Bessel equation is: 



y 



II! 



S +(i^ + 2^)IA + 1+ ( w( _^ + w _ 1) _2(li ? )) ']G(y) = , (3.57) 



V dy 



where w = Jh 3 V The solutions are: 
2(l-s) 



3-s 



(3.58) 



To obtain the expressions for the mode functions, we must identify the correspondence 
between \1/ and - as defined in (3.31) - and the solutions and This can 

be accomplished by requiring the normalizability of the vacuum. A useful asymptotic 
expansion of the two solutions is: 

4 1) ( ? /) = ^exp[%-^-^](l + 0( ? /- 1 )) , (3.59a) 

H v\v) = \[^- exp (l + O^ 1 )) , y»l, (3.596) 

Because of the exponential behaviour seen in (3.59) , we must associate ^(r],k,I) - the 

(2) 

coefficient function of the annihilation operator in (3.32) - with H„/(kr)) and ^*(r), k,I) 
with H^{krj). After we account for the normalization N and the transformation (3.48) , 
we get: 

V(ri,k,I) = N(k) (krj) w H^\krj) , I = A,B,C, (3.60a) 
ty*(rj,k,I) = N(k) (krj) w H$(krj) , I = A, B, C . (3.606) 

The alternate connection, namely \I/ <-> and <— > would lead to a vacuum 

state that is not normalizable. * 

The normalization of the mode functions comes from (3.33). When we take into account 
relations (3.60) and the Hankel function identity: 



4 2) (y> ±4%) - 4%) U?\ y) = ii , 

dy dy Try 



(3.61) 



* The normalizable vacuum state that follows from prescription (3.60) is the Bunch-Davies vacuum [4]. 
Although this works fine for power laws, there is an ambiguity for a generic scale factor. Since the solutions 
in the latter case are not known it is not clear whether a proper vacuum state exists in general. 
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equation (3.33) gives: 

s 1— 3s 

N(k) = T)Y~ S k 2 (!- s ) , (3.62) 

so that: 

V(ri,k,I) = iy/Wfj n-^T]) H%\kri) , I = A,B,C , (3.63a) 

**(77,fc,/) = J0F^ H-^ry) tf^M , I = A,B,C . (3.636) 

The basic parameter values of the most physically interesting power laws are displayed 
in Table 1. The kinds of Hankel functions generated by inflation, matter, radiation and 
flatness are presented in Table 2 and their functional form in Table 3. The explicit expres- 
sions which the mode functions take in these cases are shown in Table 4. In the case of the 
inflating universe (s = +oo) they are in agreement with previously obtained results [5]. 





Flat 


Radiation 


Matter 


Inflation 


s 





1 


2 
3 


+oo 


si(n) 


1 


V 
Vo 


V 2 
nS 


Vo 
V 



Table 1: The power law parameter s and scale factor fl(ri) for some special spacetimes. 
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V A 
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V B 


1 


3 


5 


1 
2 


v c 


3 
2 


5 
2 


7 
2 


1 
2 



Table 2: The value of the index of the Hankel functions for the special power law scale factors. 
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(9\ 

iiy (/cry) index 


(2) 

Functional Form of (fcr?) 


"-■2 


V^7 GXP 


u - 3 


V ^? exp [ ~ lk ^ 


" = i 


VX 6 ^ ^ i 1 krj (4)0 


" = 5 


6Xp [ j N (l *| (^+(^) 



Table 3: The functional form of for the needed half-integer values of the index v\ the form of 
is obtained by complex conjugation. 



Power Law 


Functional Form 


Mode Function 


s = 


vk exp [ ^ ] 


^(r),k,A) = ty(r),k,B) 


s = 


^ (l + tkrj) exp{ zkrj] 




« = * 


vk * exp [ lkr)] 




« = * 


^2= ^ ( 1 + ikrj ) exp [— ifcry] 










« = § 


2 

* ( 1 + ) ex P [ »N 






^^(i4F + t- 1 ) expHfc "i 




« = § 




*(v,k,c) 


s = +oo 


i ( 1 + %kr i ) ex p t i N 


*(ri,k,A) 


s = +oo 


V2k I 6XP [ 


*(r,,k,B) = *(ri,k,C) 



Table 4: The form of #(77, k, I) for the special power law scale factors; the form of k, I) is 

obtained by complex conjugation. 
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Step II. We proceed to calculate the propagators iAj(x;x') by substituting relations 
(3.63) in (3.43b) to obtain: 



dk exp(— ek) k sin(/cAx) 



| O(Arf) H^{kri) H${krf) + 6(-A V ) H^ikrj) ffgW) } . (3.64) 

An infrared divergence coming from the lower limit of integration may exist here. When 
present, it is due to the infinite size universe imposed by the integral approximation to 
the mode sums we have assumed. This divergence is not physical and has appeared in 
previous studies of the graviton propagator in de Sitter spacetime [6]. * On our finite 
spatial manifold the propagator is a mode sum over discrete momenta, so it is perfectly 
well defined. The integral approximation to this mode sum has lower limit kg ~ R _1 [8]. 

It is interesting to examine the range of power laws for which this infrared divergence 
exists in the basic integral (3.64). The dominant behaviour of the Hankel functions for 
small k: 

4 2) (y) — {y ~ l)Xr y~ v , (3.65) 

m 

translates to the following behaviour for the propagators: 

r+oo 

iA yi {x;x') « / dkk 2 ^-" 1 ) , fe<l, (3.66) 

Jo 

and, therefore, the following condition for the presence of an infrared divergence: 

\-vj < -\ . (3.67) 

For the "A" modes (3.67) implies: 



V A = \ 



l-3s 



s > | , (3.68) 



* For de Sitter spacetime, this does not preclude the appearance of infrared divergences in physical quantities 
[7,8]. 
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while for the "B" and "C" modes we get: 



V B 



1+8 

3-s 



\ < s < 2 , 



< 8 < I 



(3.69a) 
(3.696) 



Since the physical polarizations of the graviton are "A" modes, the corresponding crossover 
value s cr j t = | may have some indirect physical significance [9]. 

We now turn to the computation of the propagators iA l/I (x;x') for the power laws 
of interest. As Table 2 indicates, four Hankel functions are needed. Substituting the 
expressions of Table 3 in (3.64), we get for v = ^ : 

1 



iAi(x;x') 4lT 2 n{r]) n{r]/) Ax J Q 



dk sin(fcAx) exp —ik(\Ar]\ — ie) 
1 



4vr 2 Q(rj) O(V) (Ax - | Arj\ + ie) (Ax + | Arj\ - ie) ' 



(3.70) 



for v 



r+oo 
Jkn 



n( v ) o(V) Ax J ko 

iA i (x; x') - 87r 2 0(r?) " 0(r/) ~, { 



dk sin(fcAx) exp —ik(\Arj\ — ie) 



1 + 



1 + ik\Ar]\ 
k+ rj rj 



-2 + Ei 



ik§(Ax — \Arj\ + ie) 



+ Ei 



-ik Q (Ax + \Ar)\ - ie)j } , (3.71) 



for v = | 



iA 5 (x;x') 
5 



4tt 2 Q( V ) Q(rf) Ax J ko 



+0O 



dk sin(fcAx) exp —ik(\Arj\ — ie 



l + ifc|A?7| 9 + 9ik\Ar)\-3k 2 (Ar}) 2 

J- + 3 ~?j ; h 



^2 ^ ^/ 



/c 4 ?7 2 V 2 



-2i A ! (x; x') + 3iA 3 (x; x') 
1 

+ 



• » 9^ 2 - 4^(Ax) 2 + U&V) 2 + I (Ax) z - (A77) 
8tt 2 0(7^) t? 2 O(V) v' 2 1 

x I Ei[z/c (Ax- IA77I +ie) + Ei -ifc (Ax + | A?y| - ie) } j , (3.72) 
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and for v = \ 



n( v ) n( v >) Ax j ko 



P+OO 

/ 



dk sin(/cAx) exp — ik(\Ar}\ — ie) 



1 + 6 1 + ik \ Ar l\ + 45 + 45z/c|Ar/| - 15/c 2 (At7) 2 



fc 2 



7/ 77' 



/c 4 T] 2 r/ 2 



+ 



225 + 225zA;|A?7| - 90£; 2 A?7 2 - 15z£; 3 |A?7 



3 i 



k 6 r? 3 t/ 3 



hi A i (x; x) — 9iA 3 (x; x) + hi A 5 (x; x) 
2 2 2 



+ 



1 



16vr 2 0(r?) 77 3 0(77' 



n ^3 



225 /cq 



75 (Ax) 2 - 45 (At?)' 



+ 



[if (Ax) 4 - l|5 (A*) 2 (At7) 2 + If (At?) 4 ] - f ((Ax) 2 - (At?) 2 ) 

x I Ei|^o(Ax - \Arj\ + ie) + Ei -ik (Ax + \Arj\ - ie) || , (3. 



73) 



The integration technique consists of performing as many integrations by parts as 



needed to reduce the original expression to the following basic integrals: 



r+00 

Jo 



dk sin(/cAx) exp —ik(\Arj\ — ie 



Ax 



k 



^ cos(fcAx) exp 

rv 



(Ax - \At]\ + ie) (Ax + |Ar?| - ie) ' 
ik(\Arj\ - ie)] = -^j Eij^ (Ax - \ At]\ + ie) 



(3.74) 



+ Ei 



-ik (Ax+ \Ar]\ -ze)] } . (3.75) 



The first of these does not exhibit an infrared divergence but the second does and, therefore, 
its lower limit of integration is regulated to the value ko ~ R _1 . The exponential integral 
function Ei(x) is defined as [10]: 

r+00 

Ei(-x) = - dt r l exp (-t) , (3.76) 

Jx 

and for small x is very well approximated by: 

Ei(-x) ^ 7 + lnx , x<l , (3.77) 

where 7 stands for Euler's constant. It is legitimate to employ this approximation in our 

case: 



k 



-00 ,- 

— cos(/cAx) exp — ik(\Ar]\ — ie) 



k 
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^{ ln[fe§(Ax- \At]\ +ie) (Ax + |At/| - ie)] + 2j } , (3.78) 



since we have assumed throughout that koAx <C 1 and ko\Ai]\ <C 1. 
Using the approximation (3.78) and the definition: 

(x - x') 2 = (Ax - | At/| + ie) (Ax + \Arj\ - ie) , 

we can express iA L , I (x; x') in a more economical form: 

1 1 



(3.79) 



iAi(x-x') 4yr 2 0(??) 

2 r/ ?/ 



i^^3 (x^ x') 



- In 



i\2 ' 



A^(x-x') 2 -2(7-1) \ , (3.81) 



(3.80) 



iAs (x; x) = 
2 



2^ 2 V' 2 3 



8tt 2 0(77) 7/2 0(7/') t/ 2 I (x - x') 2 71 
+ 9k Q 2 - ^(Ax) 2 + |(A?7) 2 + 67/ rf + 37 \(x - x') 2 - 2r/ 7/ 



(x - x') 2 - 2?7 ?/' j \n\kl(x - x') 2 

, (3.82) 



iAj (x; x) 
2 



4t/ 3 t/ 3 

16vr 2 0(77) 7/3 0(7/') 77 /3 I (x - x') 2 
x ( In 



t(x — x') 4 — 5(x — x') 2 ?/ 77' + 12?/ 2 t/ 2 



,2/ /\2 
fe (x - X J 



+ 



+ 2 7 J + 225^ 4 - 75(Ax) 2 - 45(At/) 2 - 90?/ r\ 
^(Axf - If (Ax) 2 (At/) 2 + if (At/) 4 
55(Ax) 2 - 45(At/) 2 1 t/ 7/ + 24t/ 2 t/ /2 



(3.83) 



S^e^ 1/7. To obtain the explicit form of the scalar, pseudo-graviton and mixed propa- 
gators for the power laws of Table 1, we only need to appropriately substitute the above 
expressions in equations (3.38g) , (3.40) and (3.42) respectively. As an example we shall 
consider the case of de Sitter spacetime so that we can connect with previous work [3,8]. It 
corresponds to s — > +00 and in this limit the scalar field action becomes the cosmological 
constant term: 

; V(<p) -> (8vrG) A" 1 , (3.84) 
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ip 2 -> 



where A = 3^ . Consequently, only the pseudo-graviton propagator exists and it equals: 



pyAp a 



rpA 
\lV-Lpo 



+ iA i (x; x) 
2 



rpB 
pv 1 pa 



(x;x') = iA 3 (x;x') 

J s=+oo 2 

= iAi(x;x') | [ ; 

-^{'«[^-^ 2 ] + 2 h- 1 »} 



+ 



rpC 

pv 1 pa 



\pv-Lpa 



+ 



rpD 

pv 1 pa 



+ 



rriC 

pv ± pa 



A 



T 

\J>V->-po 



2r) rj 
8tt 2 t] 2 1 {x-x') 2 



rpA 

V-i-pa 



| In [fc§ (x-x') 2 



+ 

+ 2 



^B 



T 

(7-1)} 



+ 



rpC 

^ 1 pa 



rpA 

^ A pa 



(3.85) 



in agreement with previous results. By using the zero modes present we can absorb the 
constant term 2(7 — 1) and set /cq = [§]• However, this property is particular to the 
de Sitter spacetime and is not shared by the other power laws of interest. Finally, the 
combination of the tensor factors (3.41) that appears in (3.85) equals: 



rpA 

pv 1 pa 



+ 



B 

y.v->-po 



T 



+ 



rpC 

pv 1 pa 



Vpv Vpa 



(3.86) 



3.5 The Ghost Sector. 



The ghost Lagrangian is obtained by varying the gauge functional given by (3.14b). 
Under an infinitesimal coordinate change: 



y'» = y» + K^(y) , 



(3.87) 



the full metric transforms to: 



, dyP dy a / ,-i \ 

so that the infinitesimal variation of the pseudo-graviton field is: 



(3.88) 



pv 



' pv 



pv 



2 ^(p,v) + 2 Vpv (lnO) 7/ 



(3.89a) 
(3.896) 
(3.89c) 
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where we have decomposed Si/j^ into a term with no tfj^ and a term linear in ^^ v . There 
is a similar transformation on the full scalar field: 



<p'(x) = <p(y 1 (xfj , 



(3.90) 



and a corresponding infinitesimal variation on the fluctuating scalar field: 



5(f> = Sq4> + 5i4> , 
5o4> = —kuj <p 

The ghost Lagrangian is taken to be: 



-KU^ (p. 



(3.91a) 
(3.916) 



£ gh = -O <5F A 



(3.92) 



As usual, the infinitesimal parameter of the coordinate transformation (3.87) is the an- 
ticommuting ghost field uj and To is the antighost field. Up to a surface term the ghost 
interaction Lagrangian is: 



4^ = -^^ 

= -kQ 2 JD^ v [ % p uj% + i) vp u p ^ + Vvu/, p oo p + 2if> fiv uP (lnQ), p 
+ k (fi 2 uA \ % v u*" + \$, v u v + V w v (In Q) ,„ 
+ k 2 O 2 0'cU^^o0,p^ P , 
and completely determines the ghost-graviton and ghost-scalar interactions. 
The quadratic part of the ghost action is: 

.(2) 



(3.93a) 



(3.936) 



D A 5£ + D B o M ° o/ 



= ZJP 



P ' 



(3.94a) 
(3.946) 



leading to the following ghost propagator: 



(x, x') = iA A (x, x') rj w - iA B (x, x') ^ V P 
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(3.95) 



It is not possible to explicitly compute the ghost propagator for a generic scale factor. 
However, in the case of the power law scale factors of Table 1 we get: 



P A P 



(x; x) = iAi (x, x) 
J s=0 2 

1 1 

~ 4^2 {x - x /)2 VHP 



(3.96) 



P A P 



J s=i 



1 (x; x) = iA i (x, x) riup - iA 3 (x, x) r? M o r] p0 
2 1 



2 

m 



87T 2 ?7 2 T] 2 (x — X 1 ) 



+ 



2(7-1)) V^OVpO } , (3-97) 



2 (x;x') = iA 3 (x,x') r] fip -iA 5 (x,x') rj^o V p 
Js =3 2 2 



4 



If] T] 



- In 



kfi(x-x') 2 +2(7-1) 



87r 2 ?7 3 r\ 3 \ (x — x') 2 

§(x - x') 2 - 2t] rj' 



4 



Vpp 



In 



87T 2 T7 4 7]' 4 

+ 4 W + p - ¥( Aa 2 + l( A ^) 2 } V> • 



kl(x - x') 2 ] + 2 7 ) 



(3.98) 



,s=+oo 



(x; x ) 



iA 3 (x, x') 7/ M p - iA i (x, x) 77^0 ??p0 
5 2 



2r/ r/ 



87t 2 ?7q I (x - x') 



2 ^MP 



- In 



A: 2 



+ 



(x - x') 2 ] 
2(7 - 1) ) % p 



, (3.99) 



Again, when s = +oo , there is agreement with the previously obtained result in de Sitter 
spacetime [3]. 
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4. Corrections to the Newtonian gravitational force 

An elementary application of the perturbative tools developed above is the computation 
of the linearized response of the gravitational field to a point source. We may as well 
choose a matter dominated universe so that in the appropriate limit we can deduce the 
gravitational force law and compare it with the Newtonian law we currently measure. 

The gravitational force is obtained from the geodesic equation the worldline x M (r) of 
a particle obeys: 

r(r)+r%[ X (r)] X P(r)x' J (r) = . (4.1) 

Any deviation of its worldline from the geodesic trajectory it would occupy in the absence 
of gravitational sources will be attributed to a gravitational force. To make contact with 
the Newtonian theory we shall consider the slow motion of the particle in a weak stationary 
gravitational field. In that case, the geodesic equation simplifies considerably: 

= 204 ^00,* > ( 4 - 2 ) 
where we have chosen x°(r) = _1 (x°). The gravitational force is that part of the physical 
acceleration which is not attributable to the Hubble flow. Hence we can write: 

Fi = mn X i , (4.3) 

which, with (4.2), gives: 

F i = ^ . (4.4) 

The weak stationary gravitational field will be assumed as due to the presence of a 
point source. The action of the source worldline Xs( T ) i s: 

Ss = ~m s J dr ^-9^[x S (r)] ^(r) ^(r) , (4.5) 

and the resulting stress tensor: 

,„ s / X^r)xUr)S^(x- X s(r) 



7T(x) = -j^= / dr ^ L . (4.6) 
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The stationary nature of the source allows us to choose: 



(4.7) 



in which case the stress tensor simplifies to: 



(4.8) 



The linearized response of the gravitational field to the presence of the point source is 
given by [3]: 



jJLV 



(V,x) = ~J drf J d 3 x' n 6 (»/) [^CF*pa] (x;x) ifV) , 



(4.9) 



where the pseudo-graviton retarded Green's function is defined as: 



z^ret 
\iv^ pa 



(x;x) = 20(Ar/) Im 



pvAp a 



(x; x) 



(4.10) 



For the power law scale factors we have analyzed, the retarded propagators G T j (x; x') can 
be expressed as follows: 



GY t (x;x') = 26(A V ) Im 



iAj(x;x') ] 



(4.11) 



In terms of the mode functions: 



Gf t (x;x') = ^^- f +G ° dk ksm(kAx) lm\^(r],k,I) V*(rf,k,I) 



7T 2 Ax JO 

where we have used (3.43b). 



(4.12) 



Inspection of the form of the three tensor factors (3.41) and the stress tensor (4.8) shows 



that only the 5® 5 V Q 5 p ° 5 a ° part of 
response. 



pv-L pa 



can contribute to the desired ki/jqq(t],x) 



V>OoM = ~J drf J d 3 x' Q 6 (r/) Gf(x;x') + ^ G^(x; x')] S p ° 5 a ° ZfV) . 

(4.13) 
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Specializing to the matter dominated universe leads to: 

9(A V ) 4 f+°° 



^rct 



2 ( x 'i x ') 
s= 3 



2tt 2 Ax ^3 



/ 

JO 



dk sin(/cAx) < cos(kAr]) 



1 



A77 



(4.14a) 



+ sm(kAn) 



G 



ret 
C 



(x; x') = /" +0 ° dfc sin(fcAx) x 

V ; 2tt 2 Ax rP r/ 5 7o 



2tt 2 Ax r/ 5 r/ 5 7o 
cos (A; A?]) 



(4.146) 



225Ar? 15(Ar?) 3 45?y r/ Arj Qrj 2 n' 2 At] 



k b 



k 3 



fc 3 



k 



+ sin(kAr]) 



225 90(At?) 2 A5nn' 15n rf (A??) 2 6r/ 2 ?/ 2 3 /3 ] 



+ 



+ 



rj f) 



ju6 u4 lA £*2 £*2 

Ay Ay Ay Ay Ay 

where we have used Table 4 to substitute the appropriate functional form of the mode 
function and its complex conjugate in (4.12). We then perform all integrations by parts 

needed to express (4.14) in terms of the following integrals: 

r+oo 

is in = / dk sin(fcAz) sin^Ar?) = \ S(Ax - Arj) - %5(Ax + An) , (4.15a) 

Jo 

r+oo ^ 

I cos = — cos(AAx) sm(kAr)) = £9 (Ax + An) - Z6(Ax - Ar]) . (4.156) 

Jo k 

The result is: 



A 



/oret 



_ 2 (x;x')- ^-^{ 



S=3 Z.7I- 77 



2tt 2 „ 2 r/ 2 1 Ax ^ sm ' r? 77' 



cos 



(4.16a) 



,,2 ,/ 2 1 A.r 7 f 



6 15(x — x ) 



/\2 



?7 rf 



2r] 2 v >2 



+ 



15(x — x' 
8r/ 3 n' 3 



'COS ( 1 



which finally implies: * 

O(Arj) ijj 



/--(ret 



Gret 
C 



-i 



2 

s =3 



(x; x ) 
(x; x ) 



4tt ^2 ^/2 
0(Ar,) itf 



4tt ^2 ^2 } Ax 



— 5(Ax - An) + —7 6>(A?i - Ax) 1 
Ax nrf J 

6 15(x-x') 2 



(4.166) 



(4.17a) 



— 5(Ax - An) + 



n rj' 



+ 



15(x — x 
8r/ 3 n' 3 



2r] 2 ^2 
9(An - Ax) 



(4.176) 



These retarded Green's functions could also be obtained directly from the imaginary part of the iA 7 (x; x') 

2 

propagator (3.83). 
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When we insert (4.8) and (4.17) in equation (4.13) and use the delta function present 
in Tg (a/) to do all three spatial integrations, we get: 



_ k m s 77 J 1 



7J—X 



m 



3(x 2 -t] 2 ) 2 3(x 2 -r] 2 ) 3r/ 



4r7 3 t/ 3 Zip rf ' Arf 

where x = \\x\\ and rji is the beginning of matter domination in the universe. In terms of 
the graviton field h^ v = Q^ip^ the conformal time integration gives: 



+ 



, (4.18) 



k m s _ I 2 3r) 2 3x 3(x 2 — rf) /ri — x\ 3(x 2 — ri 
hoo(v,x) = ^^ n i "-73 "72 3 H^rJ + 



16tt 



x 4r) 6 rf 



2\2 



4?7 3 rj 2 



(4.19) 



The gravitational force is obtained from (4.4): 



,1 3 3r 3r 3r / Vtrj — r 

F(rj, r) = -G m m s { — + — — ^ + -^r; — , + ttt^ o + — — n In 



r 2 O 2 ?] 2 20 3 r/ 3 2fi 3 rj r\\ O 3 r/ 3 V Qr/i 

3r 3 1 
-2QV^}' (420) 

and has been expressed in terms of the physical distance r = O x. The first term in 
(4.20) is the Newtonian attractive inverse square law. The remaining terms - one of which 
is independent of the distance - represent the corrections. To get an estimate for their 
relative strength, we apply (1.3)-(1.4) to a matter dominated universe: 



2 1 
rj = 3 tfj *3 



V0 = 3 t , 



(4.21) 



and re-express the force as: 



F(t,r) = -Gmm s {^ + 7^ + 7^3 



7 



+ Tg r 1 3 h 3 + 



9t 3 



In 



2 - 1 



1 r 3 



13 2 1 



(4.22) 



The beginning of matter domination is thought to have occured at a physical time £i such 
that [11]: 

a{t) ~10 4 t = 10 6 ti , (4.23) 



o(*i) 
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where, of course, t is the present. Using (4.23) we obtain: 

= + l (Tf + (I) 3 + i (Lf ,410,(3-9] - W (Lf } 

(4.24) 

It is now apparent that the correction terms are negligible unless we consider objects of size 
comparable to that of the observed universe. The present time t is of the order of 10 28 cm 
while a galaxy and a cluster of galaxies have typical size r ~ 10 23 cm and r 10 25 cm 
respectively. 



5. Epilogue 

One possible explanation of the discrepancy between the predictions of the Newtonian 
theory and the data provided by the galactic rotation curves, is an appropriate modifica- 
tion of the Newtonian gravitational force law at large distances. A natural way for this 
to occur would be due to the deviations from flat space introduced by the matter domi- 
nated background. These deviations turn out to be unable to bridge the discrepancy by 
themselves or to, at least, reduce significantly the large amount of dark matter required 
for that purpose. 

In retrospect, it would have been very hard for such corrections to have the desired 
effect. For if they were strong at galactic scales to resolve the rotation curves mystery, 
they would become enormous at even larger scales something which is ruled out by the 
existence of structures much larger than galaxies: there is no gravitational confinement at 
galactic scales. The corrections would have to be of the proper strength at galactic scales 
and weak at all other scales. 

From the quantum field theoretic point of view, the Feynman rules were developed 
for a spatially flat Robertson- Walker background generated by a dynamical scalar source. 
This is a much wider class of backgrounds than what is needed for the currently observed 
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universe. Their study revealed that conformal flatness is a more powerful organizing prin- 
ciple than maximal symmetry. This fact becomes especially obvious from the identity of 
the expressions for the quadratic Lagrangian in the case of de Sitter spacetime (maximally 
symmetric) and flat Robertson- Walker spacetimes (not maximally symmetric). 

Moreover, the Feynman rules obtained make possible any perturbative calculation 
desired for power law scale factors. For instance, the existence of strong infrared quantum 
gravitational effects could be studied in detail in these backgrounds. 
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